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[::\^\\''''r/ '^'Ttie present inyehtiori generally relates to . apparatus and . 
■ methods for. prbc'es'sing inulti -component ( 3C / 4C ) 's ei smic data 
based on the use., of .recei-ver functions; more .'specif ically" \* 
, two-'dimensional * receive furictionis ' It also pertains* to' the * ' /" 
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use 



e of receiver . functions to process and interpret 

• I » • * 1 « » » 1 • • 



seismic 



■ • < 



signals;' td' derive an earth, image away .from the.'. near- siur face 
structure at . tlie; receiver, location; • ■ ' • .. . .. i 



j- 
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.BACKGROUND OF. THE INVENTION -.^V: 

■ . ■ ' • ■ . '• ■ * • ' * . ' ■■ . ' > / . •'. 

The pirocessin^ and interpretation of mult i. component. (3C/4C)V 

..'"*••• .■ * . •.* ".• ."' 

: seismic data/, acquired., directly at .the seafloor . or during, 3C . 

land: seismic surveys is compromised hy the effects , that, the *>^; 

15 ■ shallow -subsurface has on the deeper reflected, wavef ield.- 

•The* near surface is generally, associated with, low, laterally 

. varying shear wave velocities' and on land, the P-wave.,;i|f - 

• • • • 

velocity can- also be., low. These" properties often' lead, to -T^—' 

' ^^ • ■ ' • ■ ' •. ■ . ■ ' ■ • 

large/ P'r 'and S-wave ■ travel time perturbations, in th^ deeper 
20 reflected wavefield which vary from receiver to receiver. In 

• • • 

.addition, there ' are " suggestions that scattering anji 
anelastic attenuation (especially shear) impacts amplitudes 
and recorded waveforms as well . 



25 The receiver function methodology has it roots in earthquake 

. seismology, where it was. developed to investigate the 
■ ' . - ' ' * • ' ,• * . • • . 

. structure of the • crust and upper mantle using multi component 

r 'vXteleseismic v-bbdy ^ data:. . Certain: .' aspects., 'd.f '-'.; receiver 

./.functions are... described for. example by ^T. Ryberg and M. , 

Int. i41' (2000)", 1-11. Reference ' to ' * 



vll 1;'"' i'-' 



\ • ■ P f I 

■ • • '< . 



.. .•: 



. . . V* -.'. . • . . ' ' - ... ^ ^ 

t , ■ , I • t • < , • • . , • • ' 

. • • • . ■' . .* • . . . . • <.,'.••.•■' 
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While tiryirig- ' to* 'iin^ • discoiitinuxtiies ' in upper * man^ using 

receiver functions;, .1 the- Change' in the 

travel time difference between .the P-wave and the PS- 

converted ..wave, as" a' function ..of ' .the angle .'of • incidence/ 

5 could no longer -; i)e neglected i T. Ryberg and M;. ■ Weber 

proposed to calculate' a' velocity spectrum stack (VSS) based 

on the change in the traveltime difference (hereafter .called 

/ moveout) in. a single/ layer to fixid optimal- velocities for 

Stacking of .receiver functions ; " Although this approach gives 

10 good, results when the structure is relatively homogeneouis 

■ : . above the converting interface., . it lacks., a theoretical basis 

when' .applied to layered media with significant other 

discontinuities above., the converting interface, i.e. when 

ray bending is likely to produce significant deviations from 

15 the traveltime difference equation for a single layer. This 

was partly recognised and, siibsequently, ..it was suggested to 

substitute the average vertical P- and S-Wave slowness for a 

* • 

stack of layers, into the equation- for. a single layer." it' 
was also proposed to . circumvent this problem by assuming .a 

20 reference model (providing an initial vp(z) and vs(z)) and a 
corresponding one-parameter family of related models . by 
multiplying the reference model by a fraction close to one, 
and obtain the difference in traveltime as a famctibn of 
slowness by ray-tracing . through the perturbed reference 

25 model. This approach severely restricts the number of free 
. parameters. Reference to aspects of receiver function can 
also be found in the published patent applications GB 
. 2384557. and wo 02/059647.. 



* » • - ■ , .* 



30. On the other hand," normal moveout ; correction • (NMO) , velocity "'V^ 
■^-<^>V L''*-^"'^^ ^ • and ..K stacking-, of-,: ■ •.multichannel data; -has'. a . . long ■ 




• -If 





1 J • • 



. •;»* 



.•'vfc£>?.^^r.v approximatibn of 

.-^;v^.^••-■:=^M^^^'V'•v the; , squared' 



powers of. offset 




denotes' .■two-way traveltiitt^V" X;.i denotes .offset ; and* the'.';' 
' 5 •^coefficient's depend on layer* thicknesses and velocities of 



v. ... - * ... . . . . ■ , >. ...«.' 



•• • 
• 1 



'. ■ 



the medium).. Hyperbolic approximations were used, truncating '/ 
■■. ' ■ . . ' * *• . . " ' ■ . * . • • . • • * *. • ' . 

*' : "•' the infinite series, after the second term and defined, the ' 

. • rmsT' velocity . as . the : square-root of . the inverse of 

" * coefficient. .C2^-'Xt . is" aiso V known /that, the .. slope ' of .Vthe. x-.ft 

. * rf * " • ' ■ » . I ' ' .... . » ' • , , . ■ ■• • • ♦ , . " • ' ,»..!.. I ■ . • • . . - • . i • . - • » 

1 J ■, • ' ^* • ■ ' , ,..»,,...,...' 

10*' Vcurve at 'x=0 yields the inverse .bf the' squared rms-velocity 

V' ." and how'.- from . the ^ rins^velbcity ' at two ■ borisecutive depth v/v * 

levels the interval velbcity between them!', can be calculated. . ' 
' • ' . . ■ . • 

.. Taner' and ' Ko.ehler . .in: . Geophysics - 34..,, (1969) / 859-881. ;\also ,:: ■' 

- . . . . . - • . . . . .; ' • 

' •^introduced the velocity spectrum, stacking technique , using a 
15 multichannel coherence measure called semblance. The work by 

, Taner and Koehler was. generalised f or- .PS.-converted- w^ve^ by, 

... . •. ' 

' . • • • • 

' Tessmer and Behle in- Geophys . Prosp. 36 .(1988) ,. 671-688, who 
also derive a Dix-Krey type fbrmula, relating rms-velocities 

• ' ■ ; •* . ■ • * ■ ' ■ . • . • . ■ 

• for !PS-;converte.d waveg ' to- products of P- and S-wave interval 



■tt.' 

t 

f: 



20 velocities in. each layer. 



In view of the above state of art, it is an object of the 
present invention to extend and improve the use of receiver 
functions to process and interpret seismic data to derive 
25 ' images of an earth. It is further object of the invention to 

determine velocities pr velocity models from such receiver 

• ■ * ■ • * , ' .' ' . ' * •. . . . . . ....... 

functions. • - . \ 



■• f.t' 
. . ■ .11 

. v.. 



SUMMARY OF THE INVENTION 



30 



•>.". t \....-r. >...'• . . '•• . . ^» • .... ... 's 



.... , ■ _ .... I 

>.'^v-'':.^ ;;X'^r.fKThisv'» invention :• systematic; . approximations ;'..ofv^ the v;->;v,., 




• • • 



■ • ■ • 

• « • 



I'. 



J. ■', t • , . . ■ . . •, I . 1 "•■ '.LI • .• . • ' • » ■ . ..... . 



;'.:>;.•..■* ..^^/'"t" ,* ..• .•■'>.^,'^^•*•,:v 



r 



••iV- . 



' iihchanged m^^ corresijoiidirig • mode-coriveirtied signal. 



.Vlv\' : • ^ as P -waves:.: /(or *; S-waves):7and' PS-waves ' • (or ' SP~waves) , 
. converted / upon* -. iref lection or- transmission at a particular 



'.interface'", in;- a Horizontally layered medium or' a multi- 

5 ..layered earth .model.. As suchy the . relations, of the present 

' ' . • • . . . « < • • . , ... 

invention can be used to approximate the moveout of events 
• in receiver functions, which correspond to traveltime 
differences between, non-converted and converted waves. The 
receiver function can* thus be .corrected .for: moveout in- the" 
10. slowness domain or . in the time-space domain. .The moveout 

♦ • • • I " 1 * 

"* correction . overcomes limitations;^ of known , attempts, to use. 
■ receiver function to interpret an earth, which attempts were 
limited to.- shallow subsurface layers or cases of low- 

■ •.■ . , • * . » . , ... • ..J*' . .. . ' . • • V, 

,. ■ ■ . ' . ■ . • . .. _ . • ... ^ 

velocity layers* (or veirtical or near-veJrtical raypaths) . 

15 - ■ ' ■ 

In. a first aspect the traveltime difference is expanded into 
• increasing even powers of slowness, leading to a generalised 
small- slowness approximation of receiver f\inction moveout in 

! • . • 

... < ' • . - ■ 

r ' 

. the time-slowness domain. Thi^ result is more general than 
20 previous results by Ryberg and Weber (2000) since it applies 
to a stack of layers j and contains their results as the 
special case of a single layer over a half space. A pseudo 
rms-velocity is defined and this leads naturally to a Dix- 
Krey type relation for receiver fimctions, relating the 
25 • pseudo rms-velocities for two consecutive interfaces to the 
product of the P- and S-wave interval velocity between them. 



»>. • 



• ' In / a;.' second .aspect^'... .a..- related. -- .expansion of the squared - . 

• • ' , ■ " ■ ' ' • 

' . ' ■ * • " I ! . 

; traveltime difference, is made .into increasing even powers of 

^ * - . • • » I — • • / , ? '* t • - ' '/ fc- J /.*-..*■ . • ' , " * * ■ ' . '. 

30 the difference in horizontal travel distance between the P- 



• t - 



V-:- v-^:.- and .PiSr 



coriver€.iscai;.:wave' ( e . g . . dt. .=Ci+.; C2 r dx +C3 ;.dx +..:;j....wTiere dt . 

■ -T':'; •.■-;■:■ : 

r . :'v denQt'es. - ,fche..-; .•.4i.f f^rence . in'.. . traveltime. between. . . the.. . PS- ..' ■ 




•Ir 



■ • I' ^ 



. ^converted 



• - " ' --''ill 



• • • •< ♦ _ • , • 




f line t i biiv' ' i s .a 1 s o 



4 I • 



. S •) introduced J••f ormed^.by two-dimensibhai ■deconvolu€i on- of ' (x; t ) 
■ ■ -doinain' vertical component data (i:e: P-wavesV out. of tii€ 



the 

data (PS- waves) . This 



: - i.:-: corresponding " :(x; t ) • ' radial , component da 
V'.:. 2D-.receiver'. function,'.- the. -added deconvoLution . in 

' 10' the PS -wave and' the' wave as a function;- of the distanqe-l'ag- . 
■• . between , them. ' The second (truncated)- ■: escpansion , hence • can jDe ... • 
■ used ' to . NMO'^^ new type " of " receiver function. ' 

' •. '••'v"; V*'SYnthetic data,:* generated with- a-, ref lectivity . code is used.* to., • 
• V illustrated the mqveout . approximations and . .the new , type of ^ ' 
15 .receiver function. »\, 



' .■ '■ The derived series approximations (and corresponding two^l 
• term 'truncations) of - travel time', differences' between PS-^ 

converted waves and P^waves. from interfaces - in horizontally;:^ 

■ 

20 layered, media, being either a function of, slowness or a 
function of the difference in horizontal distance traveled 
.. by the two phases, enable velocity analysis and subseq[uent 
normal moveout correction of events in receiver f\inctions . 
Hie derived Dix-Krey type formulas for receiver functions 
25 relate the. .(pseudo) nos-velocities for two consecutive 
interfaces in a layered earth to the product of P- and S- 
"wave inteirval velocities between them. A new' ' kind of 
' . • '::'^''rdcei^^^ is- proposed; i calculated by. two-dimensional'- 

/ ' :. :':deconvolution of /^' (x, t) domairi.. vertical component data - out of 
30 * ' (x t) radial '-^^^-a t*t1ti r^Vi rrMre^c* information about 

• ••'r'^ ^i:> v:-^;.:7<f-;3,the tr avelt im^ 



I. 




■•„.•:•■•. . . • •• ;„,■■•;■. ri . : • ■, .■.•••.v'- • .■. , 

, ^ v^^^'----^;^'-'^^*'^^-'''-^ -''^^ >;-vv 'VV', vU..-. . 

traveled .Bm' tilao-se;;- phases . • THe.i>>-stecbnd moveokt- approximation. -vi^-;?'!' 

■ • 

.,•.•,■/.•.. •f applies to';;eyents ■iri,., this .new-, type .of receiver.- function! 

• 4 . ■ 



*• • \ 

• ■ * 



. Aspects " of the .//invention include * a' method of 

' 5 .calculating the. moveout "(or change as a- function of 
.. ' • * • ' ■ . 

■-•^srowness) ' of. -the.", traveltime . difference between a PS- 

converted wave and a P-wave, reflected/ transmitted from the* 
/same ; intierface in a layered ' . earth, " by expanding/ the 

es of increasing/ 

10 even, powers of slowness. This ' approximation can -'be made 

• increasingly , accurate ', by including mipre • terms in the 

■ 

expansion.. Coefficients of higher-order terms ' caii. be 
obtained by, regression analysis, of picked * traveltime 
/ differences iil the (t./p) domain receiver functions. '. • 
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Another . aspect of the invention is a two-term truncation of 
■the -series expansion "of the traveltime difference- between a 
PS-converted wave, and a P-*wave' reflected/transmitted from 
the . same interface and .the association of the coefficient 
multiplying the second term with a pseudo rms-velocity . This 
approximation expresses the moveout of the traveltime 
differences a in terms of a single profile or parameter of 
scalar quantities with depth rather than two profiles such 
as the P- and S-wave velocity profiles Vp(2) and Vs(z) . 

Yet another aspect of the invention is .a method of velocity 
analysis of. receiver functions calculated .by deconvolution 
o^. ..crosscorrelation ...of vertical Component . data (or . any P- 
' -wave quantity) with radial component data (or any S-wave 
■ ' guanfcity) in the' ■{t,p) domain, whicii uses the above two -term 



V* 



ft 

\'' .'^V api>roxii^ ' " multichannel, /coheiryhcy: " measure, . ' /"^^- 

•-^\^--.^'r''-'t\^}^h^^ . to;:; th"6//vel anaiysi^'-'/me.thod proposed'". W Taner • ' r:.:' 

* : ' ' '.' - '» ■' • • - • • • • ■ • . * . 'i • - , • ' . ' ' ■ . ■.. . • , ^ ' . 

■*,/:'••;.•,,:.■•/.•,;• ■ r". -.m'^ ^-^ . • • : ,• v». • • ..^V--.-'^ 





. • • • - '•«!,'''*■■ ,«?\ »*'^''* . 

J . . * ' i-V. ... » . • • 



? .V"*. ' «J' 



1. ' 



. • . ..• • 

. • I'. 

■A 



4 ; • ••• 



,y ; \. ' 'v. • . .'11 .. • 1^ .' • ». • !',. 1- • ' ' f , ' I ... •' \ • ' 'r 

• • -I* 'V '•i,' • ' • . ' 



•:" V Yet • an' 



another ' aspect -of 



f f 'I 



. ■ I *. 1' * 



this ' ' iriverition is 



a ' irietliod • of 



method 



■ • " • < — — — * I — • I. ' .," ,..',« • i ' 

' ■ ' : j^'; •> V processing ; • 'the:; " Xpseudo- ) v"eioci;ti'es . produced', by . • tlie.;,,.met 

•. • • - S : of : receiv-er ' function velocity ^ahalysis, as described above, 

■ ' oi^iaxri the ' brodiict of P-'" ' and ' ' S-wave ' propagation 



»• -I 



^' velocities in . each . layer ' of the medium/, using the Dix-Krey 

' ' ' ' . . ^ , • . ' domain receiver 



type -relations .developed f or ; the. . (t p) 

. .' . ,;' , • ...... ■ . • . . , \ ' • ' ■ ; . ». • ' . . . 

. /. ■■..functions . - v . 

: 10- • ■ • . .• . .• ■• • , ■, ■ . . 



* ^ * t 



Yet '.another astie'ct . -of vthe invention', is method . of correcting. . ;'■ : ^ 

the' events' in a (t,p) domain rieceiver function for the JNMO\ 

... \ ' . . ■■ 

of the trkveltime difference usihg^. the product ' of;, the .P- , anci;- • 
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t • - 



S-wave interval .. .velocities -obtained ; using, the method 
described above and a priory, knowledge of. either the P- or 

S-wave.:.:. velocity . or .velocity . ratio ., in y. each 

' ■ ' , , • * ■ 

theoretical e^cpressions or ray- tracing. 



» ' *** 
- -4" 



layer 



20 



Another:- aspect 
two - dimens i pna 1 
deconvolution or 
component data 



of - the invention * is. a method of calculating 
(2D) receiver functions., by 2D (stabilised) 



• It* * 



crosscorrelation of 
(or any P~wave 



(X/t) domain vertical 
quantity) with the 



corresponding (x,t) domain radial component data (or any S- 
wave quantity) . 



25 



Yet another aspect- of 



of 



the invention is a method 

' • k 

calculating, the moveout (or change" as a' function. ' of • the 
•: i difference ".^^ "-horizontal:-- ; distance..-; • traveled)- of ...-the 

• * • I , 

■; traveltime ' difference' * between- a PS-converted wave and. a. P- 

*■ '' ■ ;.' i' .' • . . ' ' ..<'•• • I . . ' • . .< •'.'.I ' • < . • 

!• •■ .. .;,s';tV.. • . V'.. '.'.'•<•". •' ■ ' ..; '/.t' ■ ; ' • ■ t. <-''■.' "• "i " . " ■.' ' ' ' *•*,• ' './,■•' 

' 30 . wave ."reflected/ transmitted from . the same interface/ " by 

• ...... . ». ."■ ," 



■• r- ... 




* 

J. i.' 
•:>•.".- 



■t 



, J> • , '. - . 



> ^*:v^ '^"-r:teSdi-stance 



t r avei;ed * 'by bo th plieLs es r Sine e 



• A" . . • 



thei:i: res\aiJB. ;,is an 



. ^ />''//;vasymp . expansion. ! (around ';dxi=0);, ; including'' more .ytexitis • in 

• • ' tne expansion does not necessarily .. improve* • the' accuracy- 
.. .. .Tiie exparisi-bh. is truncated one term ' before • the- smallest 

,t . ', • , • . , • * *■ , ' *- . . ' ■ . 

5 . which guarantees that the error* in the expansion is of the 
• • birder of the first neglected term. ■ * . 



4 • 



■•I . ■ • . » 



• ■ .( . ' • I 
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20 



" Yet another* aspect of the invention' iis • :a. method- of 



•'1 



Calculating a similar . asymptotic ' expansion - as above around a • 
10' different, expansion point (e.g." dx=2 5 Orit) .whilis maintaining • 



the possibility to; interpret the first-'-few coefficients as 
the vertical • incidence traveltime and . the . rms^velocity. 
respectively. •/ ... 



f. 



Yet another aspect of the invention is a two- term truncation 

of the. series expansion of the squared traveltime difference 

as a fxinctioxi of increasing even powers of the difference in 

horizontal distance traveled by . the two phases and the 

association of the coefficient multiplying the second term 

•* • I 

with a rms-velocity . 



Yet another aspect of the invention is a method of . velocity 
analysis of the 2D-receiver functions foanned by two- 
dimensional deconvolution using the two-term- truncation of 
25 the traveltime difference as a function of the difference in 
horizontal distance traveled and a . coherency measure,, 

' • '• " ' * ■ • ' ' . 

Similar to the velocity analysis method proposed by Taner 
.and ko.ehler .(1969 ). for .E-wave. data. " 

^ • • • • 

■ . . . . • •* . •. • • 

30 Yet another aspect of the invention is a* method- of 



■". : n'''iii:vicooBBsi:nQ tHe - velocities produfded by the method: of- 'receiver* • i..'*';^ 
.:^.'„.^ function velocity ..analysis • tp obtain the product of- P-' and 

. • • # • • ' , . , \ ' '• . * .• ■ • , ' «■.".••••'••• 




'••■'"'V- -r-.-.-.. • ..!-.v .; ■:.-s...--r-v: .. ■ 




'1 I 



... J . 



■In "i^each'-. layer • of ■''fS:Ilg?-^ed^ 
...;,v^^:::^.^^^^^^ :the Dix-Krey;,:;type..' relati-ons;* . developed foH'itlid ./;(xv t.) 

. '^dOIIlaari... two -dimensional fiinctions-. 

.-.V ' • * ' ' ' ► • . . , • '. ■ • • - • . 

*.>5.^'.v Yet.'.-, another* > aspect: . of . the • invention' -is aC method of 

/>•■. ': correcting the events in* a two-dimensional receiver^ fiinction ' 

' . *; ' . ■ '"^ V'- , , ■' 

.. for. tlie...]SIMO of the. traveltime difference using the. product 

vV' ■ of .'the P-./and' S-wave int.erval velocities ..obtained using . the ' ■ 

' abdve-^descfibed method: arid 'a priori Vloiowl edge of r..^eithe.r :the 'I''": , 



',. . . . . ) 



' 1 • 



« I ■ 



10 . P- or S-wave..' velocity 6'r velocity/ ratio . and theoretical 
expressions * or . ray-tracihg./. 



\ 4 



.* • . ■ 



• t . .. .» • , . 
• . • - ' 



Si . • ' % V ' 



" ;The ... methods .of ;.' the* . invention.,.. ' including . the..: / mdveout ;. ; 
..corrected receiver 'functions or. the velocities, derived f rdm '.-'i' 

* * ^ ' * ' ' * ' ■ • • 

15 the related Dix-Krey relations, can be applied to many known 

' s.eismic .. data ., .processing , methods, ... .including, velocity . 

. analysis, moveout correction. (JNMO and/or DM0), stacking and :;| 

' . ' • ' 

,-*.•.'■ .... 

•. pre- or post-stack migration;.* . / 



• • ft!-'- • 

"V: ?' ■ 



20 These and other aspects of the invention will be apparent 
from the following detailed description of noii- limitative 
examples and drawings . . 



BRIEF DESCRIPTION OF THE DRAWINGS 
25 FIG 1 is a schematic illustration of- the phases and raypaths 

• * • . ■ 

.'. involved in, • the . seismological receiver function 

/ . . • ♦ * . ' • ' * ' , ♦ * ' * • * . * ' ' «, • ' ■ • / * ' 

setting. A. plane P-wave is. incident on a stack of ri ' 
' horizontal . layers • from bjelow. - and mode-7 converts to*. 
Shear at. discontinuities "within, • and. at*, the base-, of . 

enoted by., solid . < lines,- S- 



30.. . .r.the. stack.;. P-waves are .denoted 



( 



' '.^ ' ' • " 




4 

■.' V , ' t .''»■ ' ' ' .I *. * 

•- ' ■l"'.Jt"i, • " "•..-.••••.,1 



" , «• • . " ,^ • . 

• ' 1 •«« 'l I ' • ■ « • 



I . 

f 

• ill 



• ■'. . •••■I 



• the ■• eictr^^r. time- ■ Th is'-:'-^i£a:^kted \ the 



traveitime rdiff erence carculatibn 



li-V.V;..,. ,. V 



. . . • •■•w 



FIG 



schematic illustration of 



the.- phases 



and 



5 



10 



15 



. FIG 



« w 



20 



2 xs a' 

raypaths . involyed in thie * reflection seismic receiver 
• function setting. A P-wave is ^ incident on a stack of 

. ' ■ • - . .• , . . , . . ... 

.n horizontal layers from above and mode-converts to 

• • • • 

Shear ujpon. reflection at discontinuities within and 

• • * . ■ : . .. ■ 

.at the; base.' . of • the • stack.. P-waves are denoted by 

■ ' t 11"' •> 

. solid -..lines ' .S^waves , by ./ dashed . • lines . The data '.xs . 

..... • • ■ _ •■■ * . 

. recorded on an .array pf equidistant multicomponent 

/receivers and the'- difference dx" in horizontal 

.♦ ^ ..... ' • 

'distance traveled by^ the P- and PS-converted *' wave 

from the same reflection/conversion point can . be* 

'••*■• ■ ' • • • V. ' ' . - • . 

. ■ " - ' '» 

'tx:eated: es:plicitly/. analogous to the difference in 
traveitime between the two phases. 



3 shows pre-processed, reflectivity data for the six- 

" ' . -• ' ■■ ■ ' ' • • 

layer over a* h^lf space model (medium properties in . 

Table 1 below) . The primairy P- and PS-converted waves 

have been identified by ray- tracing through the model " 

and are indicated in blue and red respectively. This 

is the data input for receiver function calculation 

in the slowness domain. 




■ • it' . 







9 

25 FIG. 4A are receiver functions obtained by stabilised 

deconvolution of the upgoing P-waves out of the 
. .'...upgoing PS-converted waves shown in FIG. 3. The. 

theoretical traveitime differences between the. PS- ' 

*. . * . 

converted waves arid primary P-ref lections are". 
'•."30 '^'ir!..' •.•."••» calculated* by .'ray- tracing 4rid 'shown -in blue'.'. The'. 

two-term, approximatxon (eq; 5.)..- .using exact pseudo 



. .... 



• * . . ■ 



I 




• I . • , . • • . ' . , ■ I » ■ • • ■ 1 " 



'V ' rms -velacit'ies. "1 s-'-" shown- • i'nV:Gyaii^? ' The blacfc-'tooxv-v • 
. J iV- ■■: ••■ •••^-~^^es-:tHe'-.area. wliich ,is:v'^Kbwi'';:in-''-detail . 

•• . - ; ' ' . • i f -^^ . * '. . -v. .'• ' ' ' * . • 'i. .» .. - 



•denotes 



I . ..• v."." • 



* M 'j' " 

.. •.■•. FXG. 

1 ^ / • • • " • 



r ' ....... . . 



4B are ISIMO. ' corrected receiirigrr'.f unctions - using^^^^^ two- . - 
r'term approximation (eq. 5*K-' and. the linearly;.;. . ... • 
interpolated ..exact pseudo - rms-yelpcity model. .The • 
. predicted residual mpveout is shown in blue. The 



, black box.' 



i;-^^^^. V-'^ .^^in FIG 



box. 'denotes the- area .'which, is shown ^ in detail 



10 



■ , ■ . ■ ■• , • *. ■ *- " - ; • . . • 

*/• # ,,4*1 * * 1 

FIG. 5A is a zoom- in on the area' denoted by the., black box . in. 
; . , . FIG ^A. The .two main, events correspond to the., / . * • ' ;.:f:- 
.-^ci'l at-TTra* TnrkTT-iariiii- K^^ut»^^ti ^In^i p-iref I'ect and its PS-''?v 



relative" moveout Between 'th^ P-refleqt ion and its PS 
conversion at the fifth and 'sixth interface. It is 



15 showing receiver functions obtained by stabilised • ^} 



deconvolution of the lipgoing P-waves. out of the 
upgoing. .PS-converted waves shown in. FIG. , 3 . The • '^f. 



.'I 



T - 



- •: . ■ theoretical traveltimie differences between the PS- '^5? 

I-'. I I 

i < ^ • . ^ ' ' ' , * . ' ' ' . ' • 

converted waves and primary P-ref lections are 

4 « 

20 calculated by ray- tracing and are shown in blue. The 

two-term approximation (eq.^ 5) using exact pseudo 
inms-velocities is shown in cyan. 



FIG. 5B is a zoom- in on the area denoted by the black box in 
25 FIG- 4B. The two main events correspond to the 

• ■ , ■ ■ * . * • 4 * r . « 1 ' " ' *. 

relative moveout between- the P-ref lection and its 
.-. S- .- ; {":■ ps-J'c6nv6rsioia''af"' t&e f if tfa "emii-'sixtji-'lrit^^f acei';'v It '''f^-" 



. . : . : ■" : "is showing NMO corrected^, re.cei.ver fimctioris,. using ■ ■. .;• 

• • . -. •. : "■• .' . ■ . . - . ." ■ • • • • '.' ■ '■';■>:• . '• 

the tworterm. approximation- and the . linearly .. . . • . •• • ... .. 

• -.K'. ' ihf erpbl-& t ed:«exac t pseudb: 'rms::Yeloci.ty model .v/fhe?:, ■ .- ,•: V; ■, 



', 30 



«»•»■• • 



I < 




; ;^:•;'piG•.:/ ,6A; is a receiver- functions ' calc^ by two- **'.' . . . 

dimensibnai stabilised deconvolution of (x, t) domain' 

'*/,/;•/:'- . . . '..vertical component data ojit. of .the corresponding 

. • ' ' . - . , ' ' ' ' 

5 . .radial component data. This 2D-Receiver Function" 

''''.*.'■••"■. . " '• ' . . . . " * • . • . 

• gives the traveltime difference between' a- P-wave 

* •J" * • ■ 

* • • ' » , • ^ # J,- * . , , 

reflection and its- PS-conversion as a function of the 

■ - . • . • . . ■ . ... 

• difference in the horizontal .travel distance. The 

. . ■ • . . • •■ . • - . 

theoretical traveltime differences^ .calculated by - . 
10 iraytracing, .are indicated in blue. The .two -term 

j' Approximation, calculated /f ^om eql . 20 . an^ .the .. 

* .- ' * • * • . . .• ^ • • • * 

' r. coefficients giveiivin 25a,b using -exact, medium 

properties,, is. indicated in cyan. 



i 



15 FIG. 6B shows again theoretical traveltime differences and 

the two-term' approximation are . indicated in blue, and " 
cyan respectively. Successive higher qrder 
. approximations are. shown in green. 



20 FIG. 7A is a zoom- in on the 2D-Receiver Functions and 

movebut approximations shown in FIG. 6A. The two main 
events give information about the difference in 
traveltime between the P-reflection and its PS- 
conversion at the fifth and sixth interface as a 

r 

25 fiinction of the difference in horizontal distance 

' traveled. Theoretical . traveltime differences and the 

• . ' ■ ■ . . < 

• • • . 

two- term approximation (eq.20) are indicated in blue 
and cyan respectively. 



,30,; .FIG;V 7B is. a zo6m-iri>on:. tK 



the 2D-R'eceiyer : Function^ • and 



. .... :^s:"\ \ mo veou't. approximations shown- in FIG*.* 6B. Successive '• 





■f • 

• < , . .1 

« i i 



•• v;^. ;'7^^ ' . higher :braer approximations^are also shown/ an green.- 



, -4 



■ , , • * > • \ , ■ > J . . . I ■ 'I . ■ < • , ?' , • . 
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"biETAlijED; -DESCRIPT Am).* EXAMPLES ' VV ' ' ..' . ' 



iserles expansion of the traveltime difference as a function 



■ I 



■ t 

> 



■ » 

• . ... ■■••■''««■ , . •, . 5- ' • ■ • ' • • •'*■«•.■ 1 • ■ * • •; - • • ' . -v - . 

it: >is '-fr6m.-' global'-. seisniQlogy. that, the/^dif f eaf.ence-'', 

10 ' in . traveitime . between- a P-wave incident " on ' a . stack ..of ..h • • 

•• -hori^oiitaa: layers ^ 'from- below -and. . the corresponding' ;PS-wave.; . : ' 

• • • 

converted upon transmission into the stacks can be written: 



... ■ ■ ■ - , , . ■ ■ ■ » 

1 . . 



i=l ...... * * - . . ■ 



* i ■ • < 



15: 



• ■ " : where hk ' rVsjc and Vpjc. 'denote the thickness shiear 'and| 

•compressional wave velocity of layer . k respectively/ and' 
'*•' • . • '■■ ' . • ' ' • • . . ' . ' ' • ■ ' • ' . ■ 

' the . slowness . of * the incident plane -wave i . - Note that thisfj 

formula implicitly takes intp account, the extra time (Th) it . 

20 takes the incident plane-wave to reach the horizontally • 
offset PS-conversion at the baise of the stack. Hence this 
foarmula can be used to compare the arrival time of plane P- 
and PS-converted waves with the same slowness, arriving . at 
the same receiver, as is normally done in seismology (see 

25 . FIG 1) . Equation (1). contains terms, of the form: (l-pfv^)^ ,, 
' where v can denote . either the P- or -S-wave velocity.. .Such 

• • -r-: • ' teims- can be' expanded into a- Taylor , series:- as . follows.: • 

■ • * • * 7 • • ' * . ' *' ' ' ' ' ■ 



' . • I, • ' ,>' .V. 




• .... . , • , . 



* « 



^ ■ • ■ I 



■ 



5 ' -Substituting equation ■ ( 2 ), into equation (1), the trayeitime 
difference can be. written: 



n 



• ■ , ' ■ , • - ... ... .* . . . . _ .j . 

' '"\ . . ' • " ■•.■••.*;.• \. . r — ' • . • • .• • . - . -J . ; I....- .•^/.* 

Eqpjation (4) is an expansion of the traveltime difference 
■" . between ' plane P- - .and PS- converted - waves .through . an n- layered. .• 

• > !• 

medium in increasing even powers of slowness. Note that the * 
accuracy of the expansion is directly related to the 

* 15 • accuracy, of ■ -.the ' substituted .Taylor • series - approximation 

• .« * .* 

(equation . 2) , which means that the product p^x^ should be 
small (at least smaller than 1, which means that the waves 
should be nowhere evanescent) ^throughout the stack.* 
Therefore the product of the highest P-velocity and the 

20 slowness determine the accuracy of equation 
J assximed that the slowness is small throughout 
cem truncate the infinite series given by eq. 
second term and neglect terms of fourth order jLn slowness 
and higher. This gives the small-slowness, or short-spread 

25 approximation:. 



(4) . 


If it 


is 


the 


stack. 


we 


(4) 


after 


the 



••.V 



■ It . • ' n 



• • • 



.... * . i. ^ t » • .* * « • ■ 

>• • ' . . •. •< * . ' ■ . J 



* I. V . • ' ' . J. » • . 



/ '. 30^.\.,The'. first- teorm. in edruatiph (5) .is simply,, the. difference in^ 




.1'.. ■ > 




s 



>, . ...... • ■ • . • - v. 



' V-^ ■ .difference at ' vertical incidence; dto. and the' pseudo rms- \ 

, - • *, «•' "._'»''• ',«<■ .'j ,1 ' ' . ' • ' ■ ' ■ ■ . Ill 

't'K ;s in.'Vi^^ic^jiSi^C^ap. ;respecti;v:ely,;;;fian" •.^e.;,def ii^H 'C' ' . •• rv- ;--^'/ iV-^ 

" ^ . • ' * ■ . , ■ . ■ , , , .1 • ... .1 . ■ 1 . 

10.'" • ' • * • ' 



ifc=i 

n 



V 



&b1 



• _»» 



Note that for a single layer, .the shqrtr spread approximation^ 

.given- in .equation (5) reduces to the equation given by| 

15. Ryberg . arid Weber ^ (2000). Hence, ■ their result is . a . special|. 

. ■ - • , '- * • ■ - ' 1 ■ ■ '• ■ ,1 >i*-<. 

:case of the. small-slowness " relation proposed* here. Note also:^ 

■ • • • t • * ■ 

that although equation (1) originates from earthquake 
seismology and is usually only applied in a transmission 
setting, it is also valid in a reflection setting, where 
20 plane P- and PS-waves . are considered, reflected and mode- 
converted at the same interface in a horizontally layered 
medium and recorded at the same receiver. 



5 



•■■ '-i t^^^^ between :a vei'tically' incident " P-'. and S-wa-ve (p=0); 

■ < r^^^i^^-^dvcono^: the s.t^Gkf-'bf.-'^Vlayers . Tiiy'.'.secOnd-- term/;. mult. iplyiiig. 
'V/ . is not simply ■ inteirpretable since it contains products 
'-^^^^^^^^^^^ layer ■ thicknesi^es ,and'. differences . in velocities and. has 
. /uniVs'.- . Im^/s] . , Howeyeiry • . maintaining . : th^ analogy 5tfith NMO 
' correcfci^^ developed "in • exploration and production seismic^, 

< " . ' 1 • " • I**.*.; »•*' . 

^ 1 . 1 t %' ■ ^ . . • ' . , . ■■ ■ 

■' we' c^ll ; this the . psetido ;rins-vel6city. Thus the travel time. 



•'.7 •* 



••1 : 



. ' . . < " • . ■ • . ■ • * 

A Dix-Krey relation for receiver function pseudo nas- 




V'i 



• ' 5¥'\«0S'8;8H-GB^f 3SrP. 

* . - ■ 



• ' _ . . .* 



J. 



• • • « r I" 



!• t ' I' • 4 '4 •t' >• "• ' 



• I. .• 



• 1 I •••■r.» 



. .'} 1 

.1. V. 



V n 



v. • 



-<.» 



: •Vk'^- ■siitiilar.;- ^analysi's • •a-s- 'iTessmer and'^Behlev 1988 ) 
:-;the .thiL layer k is; expressed* in 

. * bn'e=-way vertical traveltimes 



« « 



^ « 



and 



^ k 



is-.-appiied- and 
terms"' of; .the 
and ' 'iriterval 



.velocities iii that layer: 



* ■ •. 



* « 



(7) 



10 



. Also., .tke . ratio of vertical one-way. . traveltimes equals the 

' • , I i I - • 

inverse' of the ratio of in^terval velocities: ... 



« • 



k _ 



VPk 



(8) 



If equation (7) is. subsituted into the .definition - of the 
pseudo rms-velocity . (eq. , 6b) and the identity in eq. (8) is 
15 . iised, one can. derive the following . expression- for -.the pseudo 
rmsrvelocity at layers n and n-I: 



n 



n-1 



■Jfcol 



(9a) 



(9b) 



20 



Siibtracting ecjuation (9b) from ec[uation (9a) , • one arrives 
at : 



V • 



25 



~n _ ~n-l 



(10) 



• • • 1 . • • 



.using,v- thatc.the: difference ixi , the - ■ vertioal', one-way ^ 
. ... ■ S-yv/ahd P-wave . trayeltime * through. I'a^er n* is equal , to; the 

• /.'V-;; • '-.'.i^ .......... ■. - \ ; .': " 

• * ' ••< . • . • . » . • # . ,• • . . , • . I . . I ■ ■ • 



»■ 



I 




rr. . ... ^^ ... 



L I. 



' •* . . ■ 



■''■\y'6£^^^^ • s4chv;as^^^ • exp lor at: ipii^^^^ seismi-c/' 



the same, io cation 



■ the./Pr and PS-wave" energy i originating, from* 

OH' the* converting winter face are" -both* irecdrded, althbugli-' at 
. . ho'riisbnfcally offset- locations. iri\ tlie / receiver array: This 
) allows us to consider the difference in travel time, and in 

horizontal travel distance explicitly' for an array : of. 

receivers (compare FIG. 1- and FIG.. 2) ■ 



10 



15 



20 



Using elementary trigonometric, relations arid. Snell's ..law,... 
the explicit -differeriqe.' iri^'t^ a plarie^-P-r* "arid 

PS-wave • of. slowness p^ -transmitted- and converted at trie same 

... . . . . ... ... ....... : ■ ■ 

location at the- base of . layer n^ '^can be written (where hi \ 

* ■*♦.•' ■ „ . ■ 

/.VSi ^ arid Vpi denpte the . thickness, shear and compressio^al 

wave velocity of. lavei: i .respect ivelv) : " " ' • • 



k^l 



_Vs,^l-p'Vsl Vp,^l-p'Vpl_ 



(13) 



t. » 



Similarly/ the dif f erence . in horizontal distance traveled' 
can be written: 



Jt«l 



KVPk 



KVs, 



Jl-p''Vpl ^ll~p'Vsl_ 



(14) 



» > 



I- 1 . 



Now, following the analysis of Tessmer and Behle (1988) , 
both dt and dx . are e^cparided into infinite series of even 

I ' ■ 

25 powers, of p. . Hereby .the Taylor series expansion of the. 
function (l-p^v^)'^ is-, used, where v can- denote both the P- 



and S-waVe velocity. ■ 





. • '• . ^ . • N • ' 

ii • ••• t.f^ A • ^.• • . I" •♦ "^i, '« •, 



; a - J? V')~. = S]qAp^vy'-':: = l +.i.(;>^.v2) + i^X^^v^),^ +.^^^(/?^^)^ + .:..,. : (15) 



2x4 "v- ^ '• 2x4x6 ■ •■•■„'•-■.;• l^- 



a « « • I « - 

i • I 
• ■ ' ■ V 



f ' '. • . ' . p •_ !• 

I. •■ * ' 



r • 



Where the coefficients, grj gjre giyen by:" 



■ I, ^ • » i 



. • • . ■ lx3x..-x(2j-3) • 
. .2x,4x---x(27-2) 



(16) 



< ■ • 



' Note; that, this expansion is the" inverse of the expansion 

i » . « * * ,' 1 . 1 1 . . • »' * . .... . ' . * ■* - , ♦ • 

.. ' ■ ^ • ' . , • ■ ' ' 

10. used^ in , this previous ' sections (see \ equation 2) . Thus 
suibstitu ting equation (iS) into equations (13) and (14) for; 

' ^ ... . _ •'.•;.*.,■• ... . • ' ." r' 

' botH terms; containing P^'^ and , S-^ velocities to . obtain inf inite; • • ' o . . . 
series for ;dt and dx- results: in: 



n 



« « 



15 



n 



To simplify the appearance of equations (17a,b) and 
subsequent derivations, the following coefficients are 



20 defined: 



• • . >r , 

* * ^ 'I • 





n 



«».=ZW**'""'-Vpr-'). (18a) 
t-i • ■ 



in .^Zni , m+1 > • 



• I • • 



(18b) . 
(18c) - 



. .■'Usingv:thesey..v equatiorisv^;(vl7a;, b) . become : v. Vv' V . 



;fe-^-;:^>- 57:. 



r 



« 



« 



« 

I 
* 



. •■ "I 



* I 

v.. .. {* 



<^t = Xy>:P ■ (19a) 

<i*- pZ**P^ • (19b) 



>. , •• •• • 

>a)- 



• _ • < 



'•..V 



* < ' 

Equations- .19a, b have the same form as .Taner and Koehler 

•' ' ' ■ • ■ ■ • 

(1969,) .derived for the traveltime and offset of a P-wave* 
5 . reflection in a. horizontally layered medium and as Tessmer 
. and Behlev (1988). found for the traveltime and offset of PS-/ 
coriyeirtejd /wavjes / ■ suggests that . it ;is pQSB.i^ to .apply^ • ' • - 

the same methodology as they-, have done to find an expansion 
of the square, of the .traveltime difference into increasing 

10 'even powers of the difference in horizontal traveldi stance. 

*. • . . . . • . 

Such an expansion wo,uld be useful if a. new kind of receiver . 
function could be obtained/ where PS-converted waVes are 
. shifted both in time and space relative to the P-wave, 
measuring the' traveltime difference as. a function of the 
15 difference in horizontal * distance traveled by both phases. 
But thi's- new receiver function is exactly the result of 
applying* a two-dimensional (2D) deconvolution of the (x, t) 
domain vertical component data from the (x, t) domain 
horizontal component data. Thus, the objective is to express 
20 the traveltime difference dt as: 



dt^ =Ci dx^ +C3 'dx^ +C4 'dx^ , (20) 



25 The explicit form of equation (20) can be found by squaring 
* * ■ the power series for dt' (equation i9a) ,* successive higher 

■ even powers .of dx_ (equation 19b) and , ordering terms of equal 

■ 

powers in slovmess p.. Starting by squaring dx: 



■ .' , M , ' , 



• • ••• .'i^ , , ; ' - • , _ 

. • ; /' . • .> • • • ./•'k -I^-' f. •■ '.IK ; 



ft ■ ' 



M ' 




ijV; ■ 




• ^ "I 



•v: 



•..I'.. 



• • v.",- . . . 




9 I 



t-.:-.' • •'.■''''<■{/• ''.'.'.^•-f ' 'J . ... • *i jv.-i^'*.» ' • . -J ' v,-'H--v.' 

• ' ■ . • <.. . -• *<••/••'•.'. I •••• ' ■ < ■' . < " • . • • « . . •, 



' ■ *■ •••.■-V'. . * . ■ 



where , 



■ V ' f 



• • 1 



10 



> t 
i 



The 



• ■ * 



1 ' 



• '1 < 



series . • can 



be 



(22) 



found by. 



higher / ■ even power s of' thi s 

recursively applying equation' (22 ). to* find, the coefficients . , 

• • • - " ■' . -' . . , ' . ■ ' . ^ • . ' * ■ 

This leads to the . following result :' 



1 •• . .■ 



>. 'l 



2k-2 



k=l 



(23),!: ,. 



• .. . • ♦ . 

' I 



where , 



(24) 



."'1: - 



- 



Note that the coefficients Bto have to be- calculated 

* 

recursively since they contain coefficients B_,n-i- The same 
methodology can also be applied to find the square of dt to 



20 find: 



2Jk-2 



= ^1 + ^2 



K fc=i 



■ ^2fc-2 



J 



<• « -# ^ » • 



"vr 

2 



'OX ... /V. 




* 

■ . » • . •• » * 

*. • « ." .1* 



1.1, t • 
-I > . • 



where; • 



m 

- * t f 

..•VI- ' -• '. 



• • • • 



■I 



.1 



"A 



. . • '••.V' ••• 



• * • * * 

■ ■ J 



',■■<*'■ . •it' * :A>^,,'*v'',.-. . '.••.■4 ' 






• • • , « ' 



• . ■ . • . ■ ■ 

(26) , ; 



• I 



Hence, ' written out 
becomes : 



up to 6*^^ order in p, * equation (25) 



I" 



Thus the first two coefficients, using equations (18a-c),are, 



^2 = 



4j _ "^-nn _ «i _ 



5, 



11 





r 1' 






n 


r r. 






1 . . 



(27a) 



..(27b) 



Note that coefficient ci (equation 27a) can be directly 
interpreted as the square of the difference in travel time 
between the P- and PS-converted wave at vertical incidence 
(p=0 s/m) . The coefficient C2 can not be interpreted so 
easily, although the ' numerator is equal to dto, but it has 
the same form as was found by Tessmer and Behle for PS-* 
converted waves and hence it can be equaled to the inverse 
of the 'squcire of an rms-velocity . 
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AnaiogoiisV to • the derlyatibn." of . a Dix-KreY: type, ^'relation y.-f or'.* 

the*' two-term, traveltime dif f eirerice'' approximation in- the ' ' • 

* .. . *. ■ ' ' . * . . ■ * • ' 

slownesis .domain, Dix-Krey; type' relations can also be derived. 
for' .the approximations . of. the • traveltimle .difference', as a- .* 
function''- of • the' difference' in "horizontal; distance traveled - ' 



I . . I. 



by . the- P- and PS-wave. This' is the. topic of .the current •' 

section. As will be seen, the result will relate ' the rms- • 

• » ■ ' ... \. ^ ' • ■ • ■ • • ■ . 'i - .. . • • '. • ■ 

" .^••'.•V ' • ' Veloc in- the^ previoiiis section.- to" the products;/--. 

' * • . \' '"of:' P- , and;. ■ S~waye ' interval'' verocities-':' Statrting from the; '\ 

':V " ]T/.'.,Vsguared rms-velpcity : 



'.V, 



15 



■it 

*■•» 

V > 
• - .1^ - 



• aii'd . proceeding... by siibs ti tut ing /equation -'(7)'* for:'Vthe 
. thickness- of a layer - k and .the ratio .of one-way vertii|t:al 

.A* ; 

traveltimes (equation 8)- intd equation (28) results in: 

■ 



• • 25 



Hence, combining equation (2 8) and (29) , for layer n and n-1 
respectively gives : 

• « . 

. ■ 

> ■ " 

^ifrV , =y""'('rf -7f)-Vb.Vj. . .: . (30b)- 



1, , ., (.'('■" 
•• -■ • »,•'•■, . ■ . ■ 



..v.- ^.•. 
•i^?„• . 



r ' • ■ 



■ ••»..*'•»■ I • • »» ' • .« I .1.1 • » ' 1* * # ... I 

. . • ■ • • • . . ■ ; 

.' . . . . • . ■• r .. f ' , • -I ». ■ ' r »•!_. 



. . , 



r. 



■ ** ' • *n ' ;* 

* J • 



■ » 



- t t 



t . 



•l' 



NOW' finally, using that: 



5 (<-<)=;^o"-^o'' . 



(32) 



.1 •••• 



10 



eiopiatioii (31) 'can be written: 



*^ 0 mjr,n-l 



4 



: (33) 



15 



Equation* • (33 ) , is the desired Dix-Krey type . relation. • It 
shows^ how. the product of P-^ and S-waye interval velocities 
for a layer can be calculated once the rms velocities and 
the vertical incidence time differences for ■ that . layer and 
the previous layer, are- known. Equation (33) has^ exactly the 
same . form as found previoiisly for PS-converted waves by 
Tessmer and Behle (1988) . Note however that the vertical 
incidence traveltime differences and rms velocities are 
defined differently, as shown in the previous section. 



20 



Numerical examples 



25 



In the following it is illustrated how the derived 
approximations may be used with some numerical examples. 



a 



Hoveout 



two-term 



•; 1 •••• /•'• 






•'w?V^>-' . ■ 'layeiir^ a iialfspace.' The medium - p^ 

f;- suinmarisedv:;.!^ /l'.v" in- the." synthetic;/ fe>cperxirients.,.;^ 

^i''^ K-y r • ■ . \. v ^ v'i* ■ • ' ' iz^^ . ■ • * .; ■ • '■/^Jr^-- * • • 

receivers . were" positioned ;syrametrically* 'alDO-dt the sotir'ce 



- I. .■ 

• ; . .^. 

* ■ • 



••I ' 1 

■• . > 

• ■''1 , - v 

... I ' 



' position '•'With, a receiver-receiver spacing 'of ^ 3.125m, Thus,- 
5 .the fnaximxam ''source-receiver offiset cdnsidered here is 1250m.. '-../. 

■ The sampling • rate was 2ms."' The source in the modeling was a • / 

.* ■■ ' ■,..-•'...*»< ■ • . ' • . 

■ P-wave,: point, source, positioned on the surface, emitting a , 

- • . ... . ■ ■ • I . .■ . ' * 

SOHz Ricker wavelet. -.As . prd-processing step the data were . \ . 

■ • ■ . • . • ■•' . . . . , ' ' , . ' ■ " ■ . ' - ■•-*■•.'.•.,•. 

■ transformed^ to; ' ,th^ .time-^ slowness ' (xvp).- domain- The pre- ; 

10 processed: reflectivity data is shown in FIG .3. • , ^ 



r * I 



The' . ' receiver- ' ' functions, • calculated - by /stabilised 

. ■ ' . ' ' ' '"-v 

deconvolution., " for . this ....data are shown in,. FIG.;. .4A.;i*:The* 
theoretical traveltime differences between the PS-conv.^i:ted 



r 



15 waves and primiary P-ref lections are calculated by/;ray- 
. ..tracing., and .are.-, shown .in blue.. The two- term approximation 

■ ;.■,'••*' ' . .. , ' • - " ' •.>;'- 

. (eq. 5) using exact pseudo rms-velocities is 'showh in -qyan. 



It is clear that, the two-term approximation is exact * at v:^ero 
.slowness {p=0 s/m) Moreover, the accuracy decreases. j3with' 

20 both increasing- slowness and depth. This is because the 
- • • • • . . ■ . • • . 

velocity increases with depth and the accuracy of the 

expansion is directly related to the product of the velocity 

and the slowness in each layer. In FIG. 4B, the receiver 

functions have been NMO corrected using the two-teirm 

25 approximation and the linearly interpolated pseudo rms- 

* velocities. It is clear that all identified events have been 

• « * ■ ' 

. significantly flattened. • The predicted residual irioveout, 
calculated by., extracting the two- term approximation from the 
.' . ^ ray-traced; (exact) traveltime differences, is; shown in blue, 

. "' • " -■' • . •„.*.' V' . .. /. '1 I • , • ^ ■ . ♦ ' ,'r 

' • *"« 



30 



'S5^iv..A:v' ' .■In. ElGs.:'::;5^^^^^ denoted/by-;. the black::bpxes in FIGs,4s.; 



■ . . . . 



■ — ' . V. .r;-,'"'' ■ ■ • • . • . 

' • ■ < • 

. • • • 

• *• • .»#-••* 

^ ' ' ».,',•'.• *• ■*'. t'' ix , . i"''V''" *' " '•••'*''■ ' 

rela:t£vev-.. irioveout between'' the-. , P-ref lection' and its •.:-..EIS.5?:.' • 
conversion • at .* the fifth.:' and: .sixth interface. Again, it. is' . 
clear" that- the events 'have been significantly flattened 
using the two- term approximation. 



10 



i *.* •••• ■*•*, 



Note that the moveout correction using ' the two-term 

' I * ■ * » • 

approximation' obviously requires a pseudb ntis-velocity 
model. In the example here, the exact, pseudo rms-velocities ' 

is calculated . using, eq; ' .6b and; Table 1. .The. two-terin . 

.•■...*,•...■ . ^ • * . ...... 

approximation of the traveltime difference makes the problem 
amendable to velocity analysis . (Taner. and Koehler, .1969) by 
re.ducing, the number of parameters . to two. 



Layer 


P-velocity [m/s] 

» 

1 


S-velocity 

• 


[m/s] Depth [m] 


1 


1650 

J 

4 


400 


0 


2. 

• • 


1775 ■ 

■ 

■ 


700. 


■ 50 . 


3 


1900 ■ 


950 


■ 125' 


4 


2000 


« 

1000 


2 2 5 


5 


2250- 

■* 

• 


il25 


■ 425 

» 


6 


2750 


1375 


700 


7 


3000 


1500 


1300 



Table 1 Hedlum propez^ties o£ the 6-layer over a halfspace 
15 model used in the xreflectivity calculations. 



Example 2 : . Sdoveout approximations o£ the traveltime 

■ * 

difference as a function of the difference In horizontal 
20 distance traveled 



1- . . 



'1 



In this^-' example, we use -the exactly . the same. . reflectivity ';. • 

as used in example 1, with the exceptxon that the ■ ^ • 
a:''is* not' transformed to the-." time-slowness : (.T,phT''r:. 

• '• • . . , ( ;'•*... ■ .• f'i" • 

. ■ • • . . . . * . , - • ■ , * * 



data as was 




« f 



domain: v Instead, ■;• as proposed • xn • the section, on- series 
approxima'tibn: b,f ;;,the traveltiine • difference, -as; function of*; 
the difference .in horizontal .'difference, we calculate a hew' 
type of.; receiver, function '< formed . by two-dimensional (2D) 
. 5:. deconvoliitioh "of the . vertical, component (upgoing . P-waves) • 
■ out of . the radial component' (upgoing -S-waves) data in the 
(x,t) domain. - Two-dimensional . deconvolution in the (x,t) 
domain . is equivalent to. division/ of each , point in the 
; corresponding. ;.".fr^quehcy-waVenun^ .(f ,k) '. domain.. . Hence, ' 

■ 10 * "individual plane. PS- wave . components are.- shifted in .time aind 

spacd^ ; relative to - the. arrival, time . and . location of ., the P-*. 

* ' , ■ . . . , ■ ■ • 

wave - components with the- same . frequency - and ' (horizontal) 

■ . . . wavenumber..- ** 



V." 



15 The result of this' 2D-deconvolution is shown in FIG. 6A.VVThe 
. ' . theoretical * traveltime dif f erenceis. as function of-. Jcthe 
difference in horiz:ontal travel distiance, calculated by- /r;ay- 
tracing, are shown for each interface in blue.. The two-.term 

/ ' . • . . I • . . ' • • 

' . ' • . ' ■ • ' ■ . ■ . "^Sv • 

• ' . ■ "•.♦Ja , 

approximation, calculated from eq. - (20) using coefficients 

20' eq. . (27a, b) and exact medium properties, is shown in cyan. 

■* ♦ * . . , _ , 

FIG. 6B shows (in green colour) successive higher-order 
approximations. The coefficients for the higher-order terms 
have been calculated from eq. (25) , repeatedly making use of 
eq. (22), (24) and (26). Although it is not clear from FIGs. 
25 6 (or FIGs. 7)., closer inspection shows that although the 
higher-order approximations . are more accurate for small 

slownesses, they break down at a lower slowness. Thxs lias to 

', ' • ■* ' , ' . . • 

■ ' do with the asymptbtio nature- of the . series approxima.tion. 



30 Although- the 2D-receiver functions in FIG. 6A contain a- lot 



*• • ^ '1 . , . II 



6f secbiidairv ^ ^events' ■ . cajihot; be' explained- ..by., the rayK 



» in'** ' • 

• . ■ • • . »• • . . . • 



: v.„^/.,„.,.. ;.; traced V theor.etical . : traveltxme,/ di f f erences • for..,;; each . of the . ,;. 

v-r....,, . . .•.■-.■-..••■••V--;:,-.... ^K'-y' 



. . ,. • - .■'inT-ty.r • .', . 

...... . w*-'..* AVt 'Vr v 

til •Hj.tf*'-' <i'Vj 



4 • 



■ I • 1 • 

, • • . CVi," ,<f . 

: : . .■ .1 



« . • 



■^'V'>^^^^^^^ jLnt:er£aces;/ xt xs . clears.. that those events, that can. JD.e*. 

: .7; -related ; t^^^ priitiary P-ref lection ajid itis. 

• PS'-convei:sion,.' are approxiihated reasonably by the two-term' 



. 5 



apprbximatibn/ especially at small-slowness. 



In FIG. 7A, • two events from the 2D-receiver function 
corresponding to the fifth and sixth interface, are shown in 
detail. EIG. 7B shows the corresponding zoom- in on the 
theoretical/ .two-term and higher-order traveltime difference 



10 • curves . 



* • 
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20 



1/fliile the invention has been described in conjunction with 

the exemplary embodiments described, above, many equivalent 

■* ■ ■ . , , . ■<•.,'• . ... 

modifications and variations will be apparent to those 

• • • 

skilled in the art when given this disclosure. Accordingly, 

■ 

the exemplary embodiments pf the invention set forth above 
are "considered to be illustrative and hot ' limiT::ing. " ' Variou^^ 
changes to the described embodiments may be made without 
departing from the spirit and scope of the invention. 
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